Genome-wide analysis and prediction of genes involved in the biosynthesis of polysaccharides and bioactive secondary metabolites in high-temperature-tolerant wild Flammulina filiformis.

BACKGROUND: Flammulina filiformis (previously known as Asian F. velutipes) is a popular commercial edible mushroom. Many bioactive compounds with medicinal effects, such as polysaccharides and sesquiterpenoids, have been isolated and identified from F. filiformis, but their biosynthesis and regulation at the molecular level remains unclear. In this study, we sequenced the genome of the wild strain F. filiformis Liu355, predicted its biosynthetic gene clusters (BGCs) and profiled the expression of these genes in wild and cultivar strains and in different developmental stages of the wild F. filiformis strain by a comparative transcriptomic analysis. RESULTS: We found that the genome of the F. filiformis was 35.01 Mb in length and harbored 10,396 gene models. Thirteen putative terpenoid gene clusters were predicted and 12 sesquiterpene synthase genes belonging to four different groups and two type I polyketide synthase gene clusters were identified in the F. filiformis genome. The number of genes related to terpenoid biosynthesis was higher in the wild strain (119 genes) than in the cultivar strain (81 genes). Most terpenoid biosynthesis genes were upregulated in the primordium and fruiting body of the wild strain, while the polyketide synthase genes were generally upregulated in the mycelium of the wild strain. Moreover, genes encoding UDP-glucose pyrophosphorylase and UDP-glucose dehydrogenase, which are involved in polysaccharide biosynthesis, had relatively high transcript levels both in the mycelium and fruiting body of the wild F. filiformis strain. CONCLUSIONS: F. filiformis is enriched in a number of gene clusters involved in the biosynthesis of polysaccharides and terpenoid bioactive compounds and these genes usually display differential expression between wild and cultivar strains, even in different developmental stages. This study expands our knowledge of the biology of F. filiformis and provides valuable data for elucidating the regulation of secondary metabolites in this unique F. filiformis strain.

Transcriptome data reveal conserved patterns of fruiting body development and response to heat stress in the mushroom-forming fungus Flammulina filiformis.

Mushroom-forming fungi are complex multicellular organisms that form the basis of a large industry, yet, our understanding of the mechanisms of mushroom development and its responses to various stresses remains limited. The winter mushroom (Flammulina filiformis) is cultivated at a large commercial scale in East Asia and is a species with a preference for low temperatures. This study investigated fruiting body development in F. filiformis by comparing transcriptomes of 4 developmental stages, and compared the developmental genes to a 200-genome dataset to identify conserved genes involved in fruiting body development, and examined the response of heat sensitive and -resistant strains to heat stress. Our data revealed widely conserved genes involved in primordium development of F. filiformis, many of which originated before the emergence of the Agaricomycetes, indicating co-option for complex multicellularity during evolution. We also revealed several notable fruiting-specific genes, including the genes with conserved stipe-specific expression patterns and the others which related to sexual development, water absorption, basidium formation and sporulation, among others. Comparative analysis revealed that heat stress induced more genes in the heat resistant strain (M1) than in the heat sensitive one (XR). Of particular importance are the hsp70, hsp90 and fes1 genes, which may facilitate the adjustment to heat stress in the early stages of fruiting body development. These data highlighted novel genes involved in complex multicellular development in fungi and aid further studies on gene function and efforts to improve the productivity and heat tolerance in mushroom-forming fungi.

The NADPH oxidase in Volvariella volvacea and its differential expression in response to mycelial ageing and mechanical injury.

NADPH oxidases are enzymes that have been reported to generate reactive oxygen species (ROS) in animals, plants and many multicellular fungi in response to environmental stresses. Six genes of the NADPH oxidase complex components, including vvnoxa, vvnoxb, vvnoxr, vvbema, vvrac1 and vvcdc24, were identified based on the complete genomic sequence of the edible fungus Volvariella volvacea. The number of vvnoxa, vvrac1, vvbema and vvcdc24 transcripts fluctuated with ageing, and the gene expression patterns of vvnoxa, vvrac1 and vvbema were significantly positively correlated. However, the expression of vvnoxb and vvnoxr showed no significant difference during ageing. In hyphae subjected to mechanical injury stress, both O2- and H2O2 concentrations were increased. The expression of vvnoxa, vvrac1, vvbema and vvcdc24 was substantially upregulated, but vvnoxb and vvnoxr showed no response to mechanical injury stress at the transcriptional level. Additionally, the transcription of vvnoxa, vvrac1, vvbema and vvcdc24 could be repressed when the intracellular ROS were eliminated by diphenyleneiodonium (DPI) chloride and reduced glutathione (GSH) treatments. These results indicated a positive feedback loop involving NADPH oxidase and intracellular ROS, which might be the reason for the oxidative burst during injury stress.

Comparative Transcriptomics of Flammulina filiformis Suggests a High CO2 Concentration Inhibits Early Pileus Expansion by Decreasing Cell Division Control Pathways.

Carbon dioxide is commonly used as one of the significant environmental factors to control pileus expansion during mushroom cultivation. However, the pileus expansion mechanism related to CO2 is still unknown. In this study, the young fruiting bodies of a popular commercial mushroom Flammulina filiformis were cultivated under different CO2 concentrations. In comparison to the low CO2 concentration (0.05%), the pileus expansion rates were significantly lower under a high CO2 concentration (5%). Transcriptome data showed that the up-regulated genes enriched in high CO2 concentration treatments mainly associated with metabolism processes indicated that the cell metabolism processes were active under high CO2 conditions. However, the gene ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with cell division processes contained down-regulated genes at both 12 h and 36 h under a high concentration of CO2. Transcriptome and qRT-PCR analyses demonstrated that a high CO2 concentration had an adverse effect on gene expression of the ubiquitin-proteasome system and cell cycle-yeast pathway, which may decrease the cell division ability and exhibit an inhibitory effect on early pileus expansion. Our research reveals the molecular mechanism of inhibition effects on early pileus expansion by elevated CO2, which could provide a theoretical basis for a CO2 management strategy in mushroom cultivation.

The growth-promoting effects of endophytic bacteria on Phyllostachys edulis.

The research results of the growth-promoting effects of endophytic bacteria on Phyllostachys edulis indicated that the growth-promoting endophytic bacteria could improve photosynthesis in P. edulis leaves. The photosynthetic rate, transpiration rate, and the stomatal conductance in P. edulis treated with endophytic bacteria were all higher than in the control group. Endophytic bacteria could also increase the chlorophyll content and the protective enzyme activities in P. edulis, improving their reactions to the adverse environmental conditions. Through injection treatments with growth-promoting endophytic bacteria, the catalase, superoxide dismutase (SOD), peroxidase activity, soluble protein content, and soluble sugar content in P. edulis were all higher than in the control group, except for the malondialdehyde content, which was lower than in the control group.

A Jacalin-Related Lectin Regulated the Formation of Aerial Mycelium and Fruiting Body in Flammulina velutipes.

Flammulina velutipes, one of the most popular mushroom species in the world, has been recognized as a useful model system to study the biochemical and physiological aspects of the formation and elongation of fruit body. However, few reports have been published on the regulation of fruiting body formation in F. velutipes at the molecular level. In this study, a jacalin-related lectin gene from F. velutipes was characterized. The phylogenetic tree revealed that Fv-JRL1 clustered with other basidiomycete jacalin-like lectins. Moreover, the transcriptional pattern of the Fv-JRL1 gene in different developmental stages of F. velutipes implied that Fv-JRL1 could be important for formation of fruit body. Additionally, RNA interference (RNAi) and overexpression analyses provided powerful evidence that the lectin gene Fv-JRL1 from F. velutipes plays important roles in fruiting body formation.

Small GTPases and Stress Responses of vvran1 in the Straw Mushroom Volvariella volvacea.

Small GTPases play important roles in the growth, development and environmental responses of eukaryotes. Based on the genomic sequence of the straw mushroom Volvariella volvacea, 44 small GTPases were identified. A clustering analysis using human small GTPases as the references revealed that V. volvacea small GTPases can be grouped into five families: nine are in the Ras family, 10 are in the Rho family, 15 are in the Rab family, one is in the Ran family and nine are in the Arf family. The transcription of vvran1 was up-regulated upon hydrogen peroxide (H2O2) stress, and could be repressed by diphenyleneiodonium chloride (DPI), a NADPH oxidase-specific inhibitor. The number of vvran1 transcripts also increased upon cold stress. Diphenyleneiodonium chloride, but not the superoxide dismutase (SOD) inhibitor diethy dithiocarbamate (DDC), could suppress the up-regulation of vvran1 gene expression to cold stress. These results combined with the high correlations between gene expression and superoxide anion (O2(-)) generation indicated that vvran1 could be one of the candidate genes in the downstream of O2(-) mediated pathways that are generated by NADPH oxidase under low temperature and oxidative stresses.

The Sequence Characteristics and Expression Models Reveal Superoxide Dismutase Involved in Cold Response and Fruiting Body Development in Volvariella volvacea.

As the first defence for cells to counteract the toxicity of active oxygen, superoxide dismutase (SOD) plays an important role in the response of living organisms to stress and cell differentiation. One extracellular Cu-ZnSOD (ecCu-ZnSOD), and two MnSODs, were identified based on the Volvariella volvacea genome sequence. All three genes have complicated alternative splicing modes during transcription; only when the fourth intron is retained can the Vv_Cu-Znsod1 gene be translated into a protein sequence with SOD functional domains. The expression levels of the three sod genes in the pilei are higher than in the stipe. The Vv_Cu-Znsod1 and the Vv_Mnsod2 are co-expressed in different developmental stages of the fruiting body, with the highest level of expression in the pilei of the egg stage, and they show a significant, positive correlation with the efficiency of karyogamy, indicating the potential role of these two genes during karyogamy. The expression of the ecCu-Znsod and two Vv_Mnsod genes showed a significant up-regulated when treated by cold stress for one hour; however, the lack of the intracellular Cu-ZnSOD encoding gene (icCu-Znsod) and the special locus of the ecCu-Znsod gene initiation codon suggested a possible reason for the autolysis phenomenon of V. volvacea in cold conditions.

Homocitrate synthase expression and lysine content in fruiting body of different developmental stages in Flammulina velutipes.

Homocitrate synthase (EC 2.3.3.14) regulates the first step of fungal lysine biosynthesis. The gene encoding homocitrate synthase was identified in whole genomic sequencing of Flammulina velutipes and contains seven introns. The homocitrate synthase gene of F. velutipes strain W23 (Fvhcs) is 1780 bp in length and encodes a 464 amino acid protein with a predicted molecular weight 50.7 kDa. Phylogenetic analysis of Fvhcs and other homocitrate synthase proteins from diverse fungi produced a topology congruent with the current best estimate of organismal phylogeny. Analysis of protein domains by InterProScan and a motif search found that Fvhcs gene encodes homocitrate synthase protein conserved across Agaricomycotina. In addition, we sequenced the transcriptome of different developmental stages and structures of the fruiting body to analyze the expression levels of the Fvhcs gene. The data showed a correlation between Fvhcs gene expression and lysine values in different developmental stages and structures of F. velutipes.

Sequence and comparative analysis of the MIP gene in Chinese straw mushroom, Volvariella volvacea.

The mitochondrial intermediate peptidase (MIP) gene is conserved in fungi. It is linked closely with the mating-type A (mtA) gene. In this study, a fragment of the MIP gene in Volvariella volvacea (Bull. ex Fr.) Singer was first cloned by homologue-based cloning technology. Subsequently, the entire MIP DNA sequence (PYd21-MIP) was obtained after the fragment was compared with the genomic data through BLAST analysis. The PYd21-MIP sequence appeared to be homologous with the MIP gene in other fungi. Phylogenetic analysis of PYd21-MIP and other MIP sequences from diverse fungi agreed with the current organism phylogeny. Analysis of protein domains by InterProScan software and motif searching demonstrated that PYd21-MIP encodes a homologous MIP protein. These data support the hypothesis that the PYd21-MIP protein is a Hog-MIP protein homologue from V. volvacea.

Evaluation of the use of SCAR markers for screening genetic diversity of Lentinula edodes strains.

In this study, three molecular marker systems including sequence related amplified polymorphism (SRAP), random amplified polymorphic DNA (RAPD), and inter-simple sequence repeats (ISSR) were screened to select polymorphisms of 24 main commercial strains of Lentinula edodes cultivated widely in China. Twenty-nine sequence characterized amplified region (SCAR) markers were developed to set up a dendrogram using UPMGA based on nucleotide sequences of some SRAP, RAPD, and ISSR polymorphic fragments. The grouping showed that the 24 strains were apparently clustered into five groups at a level of 0.68 similarity coefficient, and those that have similar breeding background clustered preferentially into the same subgroup. Results also revealed that the 24 strains had a low level of genetic diversity, and the breeding source of L. edodes should be broadened by exploiting wild types and introducing exotic strains. In addition, the tested strains of L. edodes could be clearly distinguished and identified from others by using different combinations of SCAR primers. Thus, results of this work demonstrated that SCAR was an excellent genetic marker system to characterize and investigate genetic diversity of L. edodes. Furthermore, this provided an alternative method to identify the genetic relationship of different strains of other fungi.

Isolation of GPD promoter from Tremella fuciformis and driving expression of EGFP gene.

Glyceraldehyde-3-phosphate dehydrogenase (GPD) cDNA was cloned by RT-PCR using total RNA from Tremella fuciformis as template with a pair of degenerate primers. Then, a 500-bp 5'-upstream promoter region of the gene encoding GPD from T. fuciformis genomic DNA was isolated by thermal asymmetric interlaced PCR. The cloned promoter was fused to 5'-upstream of enhanced green fluorescent protein gene to construct T. fuciformis expression vector pCB-TEGFP with hygromycin gene as a selectable marker. Electroporation was performed to transfer plasmid DNA of pCB-TEGFP into yeast-like conidia from T. fuciformis. Molecular evidence, including PCR analysis, fluorescence detection, fluorescence spectra assay, and SDS-PAGE, indicated that the EGFP gene had been integrated into the genome of transgenic T. fuciformis strains and was expressed successfully. The results also showed that this promoter could be used to carry out regulated expression of heterologous gene products in T. fuciformis.

[Isolation and identification of a bacterial strain JS018 capable of degrading several kinds of organophosphate pesticides].

Organophosphate pesticides are used widely all over the world and play an important role in plant pest control. However these pesticides are considered as pollutants and harmful to human health. To search for microorganisms that can degrade organophosphate pesticides with high efficiency, a bacterial strain, coded as JS018, was isolated and screened from the soil in the vicinity of Shanming Pesticides Factory, Shanming, Fujian. Laboratory tests showed that the bacterium could degrade several kinds of organophosphate pesticides, such as Parathion-methyl and phoxin. The strain's degrading rates on phoxin, Parathion-methyl, hostathion and dichlorvos in LB liquid fermentation medium for 36 h were 99%, 96%, 80.4% and 69.0% respectively. The bacterial colonies on LB plate appeared shiny and pale-pink in color. The bacteria were Gram-negative coccoids, 0.5 - 0.7 microm in diameter. They grew well at 30 - 38 degrees C and pH 7.0 - 9.0. The optimal temperature and pH for cell growth was 32 degrees C and pH 7.5 - 8.0, respectively. They did not grow in medium containing 6% or more NaCl. The antibiotic susceptibility tests showed that the strain was resistant to ampicillin, penicillin and lincomycin. It was sensitive to kanamycin, tetracycline and gentamicin. Laboratory tests also showed that the strain could ferment D-glucose, trehalose, melezitose and ethanol. It was negative in the production of indole and hydrogen sulfide. It could not liquefy gelatin, utilize citrate, nor ferment L-arabinose, sucrose, D-mannitol, D-xylose, fructose, D-galactose, maltose or lactose. The catalase, urease and nitrate reduction were positive. Based on its morphological, physiological and biochemical properties as well as the 16S rDNA sequence analysis result, the strain was tentatively identified as Roseomonas sp.

Analysis of genetic diversity in Ganoderma population with a novel molecular marker SRAP.

Genetic marker technology designed to detect naturally occurring polymorphisms at the DNA level had become an invaluable and revolutionizing tool for both applied and basic studies of fungi. To eliminate the confusion on the taxonomy of Ganoderma strains, in this study, a collection of 31 accessions representative of morphotypes and some unclassified types was used for analyzing molecular diversity using a novel molecular marker sequence-related amplified polymorphism (SRAP). This collection included commercial cultivars and wild varieties that represented the great diversification of types from different countries and regions. The experimental results showed that 50 out of 95 combinations of primers turned out to be polymorphic, and 85 polymorphism bands were obtained using six combinations. Based on the appearances of markers, the genetic similarity coefficients were calculated, and genetic variations were observed (0 approximately 1) among the 31 different Ganoderma strains. The group of Ganoderma lucidum showed significant differences from the group of Ganoderma sinense. Moreover, G. lucidum in China was also different from G. lucidum in Yugoslavia. At the same time, cluster analysis successfully categorized these 31 Ganoderma strains into five groups. These results revealed the genetic diversity of Ganoderma strains and their correlation with geographic environments. It also suggested SRAP marker could be used in the taxonomic analysis of fungi. To our knowledge, this is the first application of SRAP marker on the systematics of Ganoderma strains within basidiomycetes.